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ABSTRACT 

We present observations of wavelength-dependent flux ratios for four gravi- 
tational lens systems (SDSS J1650+4251, HE 0435-1223, FBQ 0951+2635, and 
Q 0142—100) obtained with the Nordic Optical Telescope. The use of narrowband 
photometry, as well as the excellent seeing conditions during the observations, 
allows us to study their chromatic behavior. For SDSS J1650+4251, we deter- 
mine the extinction curve of the dust in the zl = 0.58 lens galaxy and find that 
the 2175 A feature is absent. In the case of HE 0435—1223, we clearly detect 
chromatic microlensing. This allows us to estimate the wavelength-dependent 
size of the accretion disk. We find an R-band disk size of = 13 ± 5 light days 
for a linear prior on r^ and r^ = 7 ±6 light days for a logarithmic prior. For a 
power-law size-wavelength scaling of oc A^, we were able to constrain the value 
of the exponent to p = 1.3 ± 0.3 for both r^ priors, which is in agreement with 
the temperature profiles of simple thin disk models (p = 4/3). 



Subject headings: accretion, accretion disks — dust, extinction — gravitational 
lensing: micro — quasars: individual (SDSS J1650+4251, HE 0435-1223, FBQ 
0951+2635, and Q 0142-100) 
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1. Introduction 

Gravitational lenses can act in many cases as cosmic telescopes. Since they can magnify 
the flux of background sources, they allow us to look far away in cosmic time, and to 
study in more detail the physics of the universe at cosmological distances, and beyond 
the resolution capabilities and detection thresholds of current telescopes. In particular, in 
gravitationally lensed quasars where variations in color are observed between their images, 
the "chromaticity" can be used to probe active galactic nucleus (AGN) and galactic structure 
using microlensing (see the review by Wambsganss 2006) and to study the properties of dust 
in intermediate- and high-redshift galaxies, as these are the effects leading to color differences 
between lensed images (see the review by Kochanek 2006). 

The extinction effects are due to differences in the amounts and properties of dust in 
the lens galaxy near each image (Nadeau et al. 1991). In many cases the extinction curve 
of the dust in the lens galaxy can be accurately determined by measuring the wavelength 
dependence of the image brightnesses. This method is similar to the "pair method" (Massa 
et al. 1983), comparing the spectral energy distributions of stars with the same spectral 
type but different reddenings, used to study dust properties in the Milky Way (MW) and 
in nearby galaxies where the stars can be individually resolved. In the MW, most of the 
observed extinction curves are well fit by the so-called Gardelli et al. (1989, hereafter GGM) 
parameterization, a relationship that only depends upon the ratio of the total to the selective 
extinction, Ry- This parameter can take values from 2.1 to 5.8 depending on the line of 
sight, with a typical value of Ry ~ 3.1 (Draine 2003). The GGM parameterization includes a 
spectral feature at 2175 A, but with some important differences that depend on environment. 
The Large Magellanic Gloud, for example, shows a weaker 2175 A feature in some regions 
(e.g., Nandy et al. 1981; Misselt et al. 1999; Gordon et al. 2003), while it is entirely absent 
for dust in the bar of the Small Magellanic Gloud (SMG, e.g., Prevot et al. 1984; Gordon et 
al. 2003). The (classical) "pair method" has been applied successfully up to distances ~ 10 
Mpc (e.g., Bianchi et al. 1996), but cannot be used at longer distances because it requires 
observations of individual stars. 

At cosmological distances, the pairs of lensed quasar images provide a powerful extension 
of the classical "pair method". Falco et al. (1999) studied extinction properties for a sample 
of 23 intermediate-redshift galaxies. They were able to estimate 37 differential extinctions, 
and also some absolute extinctions when the total absorption was large enough. The first 

measurement of an extinction curve at high redshift of comparable quality to those obtained 
in the MW was obtained by Motta et al. (2002). Mufioz et al. (2004) estimated the 
extinction laws in two intermediate-redshift galaxies and found unusual extinction curves 
compared with the MW. Mediavilla et al. (2005) extended the study of dust extinction 
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properties to the far-UV, where extinction is strongest and different behaviors have been 
observed in different environments (Gordon et al. 2003). Elfasdottir ct al. (2006) used 
theoretical analysis and simulations to study the effects of extinction in 10 intermediate- 
redshift galaxies and found no evolution of the dust properties with redshift even though, as 
they point out, a larger sample of lenses would be needed to reach a robust conclusion. Dai 
& Kochanek (2009) measured the dust-to-gas ratio in six distant lens galaxies using Chandra 
observations to measure the difference in the gas column for comparison to the extinction 
differences, and they found a value consistent with the average ratio in the MW. 

There arc few other methods for studying dust extinction outside the MW and nearby 
galaxies. Another "non-classical" pair method compares reddened and unreddened photo- 
metrically similar Type la supernovae(SNe la; Riess et al. 1996; Perlmutter et al. 1997). 
Quasars with foreground absorbers have been used to measure the excess reddening pro- 
duced by these systems. For instance, Murphy & Liske (2004) studied the properties of 
dust associated with damped Lya systems. Galactic environments were also traced with 
absorbers such as CIV (Richards et al. 2001), Call (e.g.. Wild et al. 2006), and Mgll (e.g., 
Menard et al. 2008). Finally, the spectra of the afterglows of gamma-ray bursts (GRBs) 
have also been used to determine extinction curves of dust in the host galaxies (e.g., Jakob- 
sson et al. 2004). These methods have found that extinction curves similar to those in our 
neighbor galaxies are also found at high redshift. The 2175 A feature has been confirmed 
in several systems (Motta et al. 2002; Junkkarinen et al. 2004; Wang et al. 2004; EUison et 
al. 2006), and its highest redshift {z — 2.45) detection to date was made by Elfasdottir et 
al. (2009) using GRBs. The SMC extinction law also appears to reproduce the extinction 
in some intermediate- and high-redshift systems, such as the line of sights of many quasars 
(e.g., Hopkins et al. 2004), and to other GRB afterglows (e.g., Kann et al. 2006), although 
the sample of robustly confirmed SMC-like extinction laws is very small. 

Besides extinction, microlensing-induced variability by stars and compact objects in the 
lens galaxy (see the review by Wambsganss 2006 and references therein) can also produce 
chromatic variations between the quasar images. The magnification produced by microlens- 
ing depends strongly upon the source size, with smaller emission regions showing larger 
variability amplitudes. In particular, since quasar accretion disk models (Blaes 2004, and 
references therein) all assume that the thermal emitting region is larger at longer wave- 
lengths, we should observe different microlensing magnifications at different wavelengths 
(i.e., chromatic microlensing). Although non-negligible, this chromatic variation is more 
difficult to observe compared with achromatic microlensing, because microlensing magnifica- 
tion gradients on the scale of the source size are required to produce the effect (Wambsganss 
& Paczynski 1991). However microlensing-induced chromaticity has been detected and an- 
alyzed in many lens systems (Anguita et al. 2008; Bate et al. 2008; Eigenbrod et al. 2008; 
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Poindexter et al. 2008; Floyd et al. 2009; Mosquera et al. 2009; Blackburne et al. 2011). 
The sizes of the non-thermal emission regions have also been studied using the differences 
between optical and X-ray flux ratios (Pooley et al. 2007; Morgan et al. 2008; Chartas et 
al. 2009; Dai et al. 2009; Blackburne et al. 2011). These results showed that microlensing 
is a unique tool to zoom in on AGN and measure their physical properties. Unfortunately 
it is not always easy to separate extinction from chromatic microlensing when variations in 
color are observed. 

One successful approach to studying chromatic microlensing is the use of narrowband 
photometry (Mosquera et al. 2009). Since the broad line emitting regions of quasars are 
generally too large to be microlensed (e.g., Abajas et al. 2002; Bentz et al. 2009), offsets 
should be observed between the flux ratios of adjacent continuum and emission line wave- 
lengths when microlensing is present. These offsets cannot be accurately measured using 
broadband photometry, since each filter typically contains contributions from both emission 
regions. Narrowband filters can separate the two, and thus distinguish microlensing from 
extinction. 

In this work, we examine the wavelength-dependent flux ratios of SDSS J1650+4251, 
HE 0435-1223, FBQ 0951+2635, and Q 0142-100 measured at eight different wavelengths 
and at several epochs using the Nordic Optical Telescope (NOT). These systems were selected 
because they have previously shown evidence for extinction or achromatic microlensing. The 
details of the observations and data reduction techniques are described in Section 2. The 
modeling and analysis of each individual system appears in Section 3. The conclusions and 
the main results are summarized in Section 4. 

2. Observations and analysis 

We observed the gravitational lens systems SDSS J1650+4251, HE 0435-1223, FBQ 0951+2635, 
and Q 0142—100 with the 2.56 m NOT located at Roque de Los Muchachos, La Palma, Spain. 
The images were taken with the 2048 x 2048 ALSFOC detector, which has a spatial scale of 
0.188 arcsec/pixel. We used seven narrowband filters, plus Bessel-I, which cover the wave- 
length interval 3510-8130 A. In the case of SDSS J1650+4251, a second round of monitoring 
observations was performed using the Liverpool Telescope (LT) with the 2048 x 2048 RAT- 
Cam CCD, the scale of which is 0.135 arcsec/pixel. We used the Sloan Digital Sky Survey 
(SDSS) ugriz filters covering the wavelength range from 3000 A to 10000 A. Tables [U and 
|2] provide a log of our observations and describe the filters. The data were reduced using 
standard IRAF procedures. After overscan and bias correction, a master fiat-field image was 
created for each filter and then applied. Finally, frames of the same object were aligned and 
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combined for each filter. Cosmic rays were removed when combining the frames by using a 
3cr crreject rejection algorithm. 

We used imfitfits (McLeod et al. 1998; Lehar et al. 2000) to model the images and 
derive the magnitude differences between the quasar images as a function of wavelength. The 
quasar images were fit as point sources and the galaxy luminosity profile was described either 
with a de Vaucouleurs model or an exponential disk. The model consisted of the positions 
and intensities of all the components, as well as the parameters describing the lens galaxy 
morphology (effective/scale radius, ellipticity, and position angle). As these parameters 
were allowed to vary, the models were convolved with several point-spread functions (PSFs; 
typically three or four) defined by stars near each lens system and optimized by computing 
the to fit the observed image. If Hubble Space Telescope (HST) data were available, the 
component positions and the structure of the lens galaxy were fixed to the values derived 
from the HST data. 

Since the main contribution to the residuals is caused by systematic errors arising from 
the PSF models rather than statistical fluctuations, we estimated the final magnitude differ- 
ences from the model whose PSF best fits the combined frame. The final uncertainties in the 
magnitude differences were calculated by adding in quadrature two contributions: one that 
comes from the dispersion of the results from the PSFs and the other from the dispersion of 
the fits to the individual frames. In all cases, the residuals were negligible compared to the 
image fluxes. Tables E] and H] report the results, and the details of the individual models are 
discussed in the next section. 

3. The Individual Systems 

3.1. SDSS J1650+4251 

SDSS J1650+4251 is a two image lens discovered by Morgan et al. (2003). It has a 
source redshift of Zs = 1.547. The separation of the two QSO components is 1'.'2 with B-band 
magnitudes of 17.8 and 20.0. The detection of Mgll and Fell absorption lines in the spectra 
of the lensed quasar suggests a lens redshift of zl = 0.58 (Morgan et al. 2003), although a 
resolved spectrum of the lensing galaxy is needed to confirm this value. The time delay for 
the system is estimated to be ~ 1 month (Morgan et al. 2003; Vuissoz et al. 2007). 

The first series of observations of SDSS J1650+4251 were obtained with the NOT tele- 
scope using eight mostly narrowband filters. Due to the good seeing conditions, the two 
quasar components A and B (the brightest and the faintest, respectively) were clearly re- 
solved. An unassociated galaxy 4 arcsec south of image A was also detected. To estimate the 
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I-band magnitudes, we fit the two components as point sources and modeled the lens galaxy 
and the nearby "external" galaxy as either a de Vaucouleurs model or an exponential disk. 
Since the lens galaxy is very faint, it is difficult to distinguish between these two luminosity 
profiles. The best fit was obtained with a de Vaucouleurs profile for the lens galaxy and 
an exponential disk for the nearby galaxy. The excellent seeing conditions (0'.'6) allow us to 
determine the relative positions and brightnesses of the quasar images, as well as the position 
of the lens galaxy and its flux. The results for the relative component positions in Table [5] 
are in agreement with the ones obtained by Morgan et al. (2003). The broad I-band filter 
is the only one in which the lens galaxy is well detected, appearing as a clear residual if we 
use a model excluding it. 

In the seven narrowband filters, the contribution from the lens galaxy is very small, 
and the lens system can be modeled simply as two point sources. The resulting fits are 
excellent, as illustrated by the residuals for the Stromgren-y image shown in Figured] The 
peak residuals for the images are less than 1.7% of the peak intensity of image A. The 
absence of the lens galaxy in the residuals, even at wavelengths beyond the 4000 A break, 
is consistent with the expected level of contamination. Assuming that the calculated flux 
ratios (with and without galaxy) are well determined in the I-band, we can estimate the 
contamination in the bluer filters using a quasar spectrum template for the images and an 
early-type spectrum template for the lens galaxy (Assef et al. 2010). We find that the lens 
galaxy contamination in the flux ratios would be at most ~ 0.04 mag in the bluer filters, 
within the level of the model uncertainties. The photometric results appear in Table S] and 
are plotted as a function of wavelength in Fig. [2] (black squares). In Figure [2] we have also 
included the results of Morgan et al. (2003) (empty triangles), which are in good agreement 
with the ones we obtained. 

To model the effects of differential extinction we fit the magnitude differences as a 
function of wavelength following Falco et al. (1999) as 



where AM = Mb-Ma is the relative magnification, /\E{B -V) = Eb{B-V)- Ea{B - V) 
is the differential extinction, and -R(A) is the mean extinction law. The fits can also be done 
as a function of the "dust redshift" (Jean & Surdej 1998), either to estimate the lens redshift 
or as an added confirmation of extinction. All these quantities can be determined without 
needing to know the intrinsic spectrum of the quasar, mo (A). We use a statistic for the 
fits and either the CCM (1989) parameterized models for the Galactic extinction curve or 
the Fitzpatrick & Massa (1990) model with its parameters set to the values found by Gordon 
et al. (2003) for the average extinction in the SMC. The results are shown in Figure [2J 




(1) 
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Our best fit for tlie data witli Zdust = 0.58 is obtained for an SMC extinction law 
with Xdof = O.odll, where xL/ is psr degree of freedom. The parameters for this fit are 
AM = 1.8 ±0.1 and AE = 0.10 ±0.02 at la. For the CCM extinction model with Ry = 3.1, 
we find Xdof = 2.2, which is a significantly worse fit because the Str-u point is in conflict 
with the presence of a 2175 A feature. Excluding the Str-u point, the CCM model would 
also fit well (xL/ = 0.03). 

The alternative to extinction as an explanation for the wavelength-dependent flux ratios 
is chromatic microlensing. We view this as a less likely explanation because the wavelength 
dependence observed in the NOT data is consistent with that observed ~ 4 months earlier 
by Morgan et al. (2003) and ~ 3 years later in our LT observations, as shown in Figure |3l 
There is a wavelength-independent shift of ~ 0.2 mag that could be due to microlensing, but 
the chromatic structure is unchanged. Since the observations made at the LT cover a period 
larger than the expected time delay and no significant variations were observed between 
the different nights, the magnitude shifts are probably not due to intrinsic variability in the 
quasar modulated by the time delay. 



3.2. HE 0435-1223 

HE 0435—1223 was discovered by Wisotzki et al. (2000) in the Hamburg/ESO survey 
for bright QSOs and later identified as a gravitational lens (Wisotzki et al. 2002). In this 
quadruple system, a background quasar with a redshift of zs = 1.689 is gravitationally lensed 
by an early-type galaxy at = 0.4541 (Morgan et al. 2005). Integral-field spectrophotom- 
etry (Wisotzki et al. 2003) showed no evidence for differential extinction, while there was 
evidence for microlensing. Kochanek et al. (2006) also observed microlensing variations 
~ 0.1 mag yr~^. Due to the symmetric distribution of the images around the lensing galaxy, 
the time delays of At ab = -8.00^0:82 days, AtAc = -2.10lo:7i days and At ad = -14.37^0:85 
days (Kochanek et al. 2006) are relatively small. Courbin et al. (2010) found similar time 
delays using a longer monitoring period but with larger formal uncertainties. 

Observations of this system were made using the narrowband filter set at the NOT on 
two nights separated 15 days (see Tabled]), and the photometric results for the two epochs 
agree at the l-cr level. Figure H] shows the results of the PSF photometry for one of those 
nights, in which the lens galaxy was modeled with a de Vaucouleurs profile, and the quasar 
images were included as point sources. Their relative positions and the structure of the lens 



^We neglect the somewhat offset 1-band point. We think that it is discrepant because the lens galaxy was 
slightly oversubtracted. 
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galaxy were fixed to the HST values (Kochanek et al. 2006). 

The behavior of the flux differences in Figure H] looks very similar to what we found 
when studying Q 2237+0305 (Mosquera et al. 2009), suggesting that we have detected 
chromatic microlensing in another quadruple system. Considering the quasar redshift, the 
Stromgren-u and Stromgren-v filters are the only ones affected by the broad emission lines 
of the quasar. The Lya emission line contributes about 38% of the flux in the Stromgren-u 
filter, and part of the CIV emission line extends over the full width of the Stromgren-v filter 
based on the SDSS composite quasar spectrum (Vanden Berk et al. 2001). Therefore, in 
a system affected by microlensing, these two filters will be offset from the continuum flux 
ratios at those wavelengths, because the broad line regions are much larger than the length 
scales for microlensing (e.g., Abajas et al. 2002; Bentz et al. 2009). 

Looking in Figure S] at the filters that are not affected by emission lines, it appears that 
chromatic microlensing is affecting image A. The observed chromaticity between the bluest 
(Str-b) and the reddest (I-band) filters is (Am)/.;, = 0.20 ± 0.09. As none of B, C, or D 
shows significant signs of a wavelength dependence, even at wavelengths contaminated by 
emission lines, achromatic microlensing must be weak for these images. This is supported by 
the R-band light curves obtained with the SMARTS 1.3 telescope located at the Cerro Tololo 
Inter- American Observatory, in Chile (Blackburne & Kochanek 2010). Figure |5] shows the 
brightness fluctuations observed in the difference [rriB — m^), once a time delay correction of 
—8.0 days was applied, and time-dependent changes are not seen in the other two differences 
{niB — mc ~ 0.05 and — rriD ~ —0.25 ). Thus, image A seems to have been undergoing a 
microlensing event near the time of our NOT observations, while no-microlensing variations 
were observed in B, C, and D. Our flux ratios for the redshift zero Ha filter agree well with 
the similar- wavelength SMARTS R-band observation (Figure |5l filled squares). 

We modeled the microlensing of image A as follows. First, we fit a simple singular 
isothermal ellipsoid plus external shear model to the HST positions from Kochanek et al. 
(2006) using the lensmodel package (Keeton 2001) to determine a convergence and shear for 
image A of k = 0.43 and 7 = 0.39. Second, we generated microlensing magnification patterns 
using the inverse polygon mapping method of Mediavilla et al. (2006). We generated patterns 
with stellar mass fractions of = 0.01, 0.05, 0.1, 0.15, 0.2, 0.25. We used stellar masses 
of M=l Mq, an outer scale of 20 and a dimension of 2048 x 2048 pixels to get a resolution 
of 0.01 r£;/pixel. With the present data we cannot determine M, and our size estimates can 
be rescaled as (M/M©)^/^. Finally, we estimated the R-band microlensing magnification 
of image A assuming that the H-band fluxes from Kochanek et al. (2006) represented 
the true flux ratios. The H-band data are the reddest ones reported in the literature for 
HE 0435—1223, and will be less affected by extinction and chromatic microlensing. Even 
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though chromatic effects could be contaminating the H-band measurements, our observations 
support the hypothesis that they are neghgible. The H-band magnitude differences, (m^ — 
mc)H = -0.02 ± 0.04 and (m^ - m^)// = -0.23 ± 0.06, match the respective SMARTS 
R-band differences (see Figure E]). Since from the NOT chromatic data we know that neither 
B, C, nor D is significantly microlensed, the H-band flux ratios should be a good estimate 
of the real fluxes of the images. This means that at R-band image A is magnified by 
microlensing by AM^ = —0.19 ±0.04. The differential microlensing of image A between the 
Str-b and R bands was simply estimated from the observed chromatic microlensing, and is 
AMb-R = -0.13 ±0.08. 

Finally, we convolved the magnification patterns with Gaussian source profiles (oc exp(— r/2rs)) 
of varying size to model the R-band and = r^{Xiy/ X^y for the Str-b band to compute 
the probability of reproducing AMr and AM^-r as a function of the size rf, the power- 
law index and the stellar fraction /t*//t (see Mosquera et al. (2009) for details on the 
probability distribution calculations). We adopted p = ^, i =1, 6, and rf = (1 ± 2i) 
light-days, i =0, 12. Figure |6] (solid lines) shows the probability distribution P{rf,p) for 
the different values of the stellar fraction a = The contours correspond to 15%, 47%, 

68%, and 90% confidence intervals. The integrated values of the probabilities, Pa{r^) and 
Pa{p), are shown in Figures [7| and [8], respectively, for the considered stellar fractions. If we 
adopt the likelihood maximum as our estimator for and p, they can be constrained at 68% 
confidence for a < 0.15 (Table [6]). For a > 0.15 the shape of the probability distributions 
does not allow us to determine the corresponding uncertainties. However, we know from 
Mediavilla et al. (2009) and from time delay measurements (Kochanek et al. 2006) that 
low values of k^:/k, < 0.2 are favored. Therefore we calculated P{r^) and P{p) using the 
probability distribution of the stellar mass fraction found by Mediavilla et al. (2009) as a 
prior (solid line in Figures [91 and ITOl) . interpolating over the probability distributions since 
the sampling in a is not uniform. From these weighted probability distributions we estimate 
that the disk size and the power-law slope are = 13±5 light days and p = 1.3 ±0.3 at 68% 
confidence. We also estimated these values for a logarithmic prior on (Table [6]) just by 
dividing by the distributions obtained with a uniform prior. The probability distributions 
P{r^,p) using this prior, and the corresponding integrated values Pa{r^) and Pa{p), are also 
shown in Figures El U\ and [HI respectively (dashed lines). In this case we found a disk size 
of = 7 ± 6 light days and a similar slope p. Pooley et al. (2007, 2009) and Blackburne 
et al. (2011) found similar disk sizes to these estimates by comparing optical and X-ray flux 
ratios, and similar results were found by Morgan et al. (2008, 2010) modeling the R-band 
light curves (they estimated = 15^9^ light days for a face-on quasar). The slope p for 
the wavelength dependence is in agreement with the Shakura & Sunyaev (1973) simple disk 
model {p = 4/3). Blackburne et al. (2011) also found solutions compatible with our results 
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for HE 0435—1223 with p = 0.55 ±0.49 and p = 0.67 ±0.55 for a linear and for a logarithmic 
prior, respectively, although they found a shallower average slope for their full sample of 
lenses. 

If we compare the probability distribution P{a) without the Mediavilla et al. (2009) 
prior, it increases monotonically with a (Figure [TTl dashed line). This is the expected 
behavior for P{a), since a single epoch chromatic microlensing detection likely introduces 
a strong bias toward large a values. With a large there are many more regions with the 
strong micro-magnification gradients needed to produce the observed chromaticity. This 
is the reason why we cannot constrain the stellar fraction with our procedure. When we 
include the Mediavilla et al. (2009) prior, the resulting P{a) resembles that prior (Figure 
[TTl solid line). In the end, our procedure is essentially equivalent to computing the probability 
distributions P{r^) and P{p) for a fixed value of a ~ 0.1 that is also consistent with the 
time delays. 



3.3. FBQ 0951+2635 

The gravitational lens system FBQ 0951±2635 was discovered by Schechter et al. (1998). 
It has two quasar images separated by I'.'l, and the time delay between them is estimated to 
be approximately two weeks (e.g., Jakobsson et al. 2005). The quasar redshift is zs = 1.246. 
However the redshift of the lens galaxy was much more difficult to measure, since it was very 
difficult to disentangle its flux from the quasar images. Kochanek et al. (2000) suggested 
a value of zl = 0.21 from the position of the lens in the fundamental plane, and this was 
spectroscopically conflrmed to be zl = 0.260 ±0.002 by Eigenbrod et al. (2007). Indications 
of microlensing in the system were found by several authors (e.g., Schechter et al. 1998; 
Jakobsson et al. 2005), but the chromatic behavior of FBQ 0951±2635 is still not well 
understood. This system was observed during one night at the NOT (see Table [1]) using our 
narrowband fllter set. We flt the data using point sources for the quasar images, and the 
lens galaxy was included in the I-band model as a de Vaucouleurs proflle. Figure [T2] shows 
the result of our photometry compared with HST measurements at two different epochs 
(CASTLES and J. A. Munoz et al. 2011, in preparation). 

If we focus on the NOT data in Figure [12] (open squares), our results are compatible 
with no chromaticity. However, our data shed little light on the chromatic behavior of 
FBQ 0951±2635, because almost all the fllters contain strong emission lines. The iac#29 
filter corresponds to the OIII line, iac#28 is blended with Mgll, the [NelV], Felll and [OH] 
lines lie in the Str-y filter, Str-v is affected at ~ 50 % by CHI] , and the Str-u filter contains 
SilV, OIV], CIV, and Hell. In any case, our results as well as those obtained by J. A. Munoz 
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et al. (2011, in preparation; Figure [T2|) . confirm the lack of extinction because no cliromatic 
fluctuations are observed. The J. A. Mufioz et al. (2011, in preparation) flux ratios are 
shifted by ~ 0.1 mag, while the earlier CASTLES observations show a larger shift and a 
significant wavelength dependence, suggesting that there is significant microlensing in this 
system. In any case further observations are needed to better understand the chromatic 
behavior in this system. 

3.4. Q 0142-100 

The gravitational lens Q 0142—100 (UM 673) was discovered by Surdej et al. (1987). 
It is a doubly imaged quasar with components separated by 2'.'2. The source redshift is 
zs = 2.72 (MacAlpine & Feldman 1982), and absorption lines detected in the quasar spectra 
suggest a lensing galaxy at zl = 0.49 (e.g., Surdej et al. 1987). Several authors have 
discussed the chromatic behavior of this lens system (Falco et al. 1999; Wisotzki et al. 
2004; Elfasdottir et al. 2006), but its nature is still a matter of discussion. The photometric 
data obtained by Nakos et al. (2005) possibly detected chromatic microlensing and renewed 
interest in this lens system. The narrowband observations of Q 0142—100 were obtained at 
the NOT on two nights separated by 15 days (see Table [1]) . The best fit to the images was 
found using a de Vaucouleurs profile for the lensing galaxy, and point sources for the quasar 
images using the astrometry and structural models from Lehar et al. (2000). 

Unfortunately, these observations shed little light on the chromatic behavior of Q 0142—100 
(open squares in Figure [T3|) . Because most of the NOT filters are again contaminated in 
varying amounts by the emission lines of the quasar. The most affected filters are the iac#29 
and Sty-b bands, which lie on the GUI] and Lya emission lines. However, comparing our 
results with the data obtained by Elfasdottir et al. (2006) (filled squares in Figure [13]) and 
by CASTLES (Falco et al. 1999) (filled triangles), somewhat constrains the origin of the 
observed chromaticity. All three data sets, spanning a 10 year period, show the same color 
trends. Koptelova et al. (2010) also monitored this system during a 2 year period in V, 
R, and / bands, and their average niB — rriA values for the different filters are consistent 
with our observations. This similarity between the color trends essentially eliminates the 
possibility that intrinsic source variability modulated by the time delay could produce the 
color trends, since the delay, while unknown, is small compared with the 10 year interval 
between epochs. 

Extinction alone is also ruled out because the wavelengths of the observations corre- 
spond to the regime where all the extinction laws are roughly proportional to A, with the 
2175 A feature lying blueward of all the data, leaving nothing to create the observed parabolic 
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wavelength dependence. Therefore the explanation for the observed chromaticity should be 
chromatic microlensing, even though the chromaticity changed little over a decade. Mi- 
crolensing variations for Q 0142—100 are assured on the timescale — 24 year^ it should 
take the source to cross an Einstein radius, although fluctuations can occur on the shorter 
timescale of ~ 3 years it should take to cross the optical source siz^. This suggests 
that the images must be lying in the broad, relatively structureless "valleys" between the 
"active ridges" of the caustic networks. The redder color of image B may partly be due to 
contamination by the lens galaxy, since their separation is only 0.38 arcsec and our data are 
limited by the resolution of our ground-based observations. 



4. Summary and Conclusions 

In this work, we studied the chromatic behavior of four lens systems using optical multi- 
wavelength data from a monitoring campaign performed at the NOT. The use of narrowband 
photometry under excellent seeing conditions has proved to be a powerful tool to disentangle 
chromatic microlensing from extinction, since depending on the source redshift, emission line 
contamination can be better separated from the continuum emission. 

In particular, the color variations in SDSS J1650-I-4251 are probably dominated by 
extinction and require an extinction law without a 2175 A feature. This indicates that we 
have found a galaxy in which the extinction curve is similar to that of the SMC but at a 
cosmological redshift {zl = 0.58). This result significantly increases the small sample of 
high-redshift galaxies with similar extinction features, and it is the first robustly determined 
from the pair method. This was possible due to the good wavelength resolution achieved 
by the use of narrow filters. The firm confirmation of an SMC-like dust at a cosmological 
distance is very important since many models including dust at higher redshift s assume the 
"featureless" SMC extinction law (e.g., Richards et al. 2003; Hopkins et al. 2004). The 
selection criteria would certainly be crucial in many fields of astrophysics, like, for instance, 
in the correct understanding of the expansion of the universe through SNe la (e.g., Jha et 
al. 2006, 2007). 



^Tliis timescale was calculated using a concordance cosmology {ilm = 0.3, 57a = 0.7, Hq = 72 km s^^ 
Mpc^-'^). The Einstein radius is te — 2.7 x 10^^ cm for 0.3 Mq microlenses, and the effective transverse 
velocity of the source (e.g., Kayser et al. 1986) is w « 360 km s~^ for a velocity dispersion of cr* = 237 km 
s-i. 

•^We estimated the V-band disk size assuming that the disk emits like a thin disk (Shakura & Sunyaev 
1973) and used the black hole mass of Mbh = 2.26 x IO^Mq estimated by Peng et al. (2006). For these 
assumptions the disk size in the V-band (Arcst — 0.15/xm) is Ry — 2.9 x 10^^ cm. 
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We clearly detect chromatic microlensing in the image A of HE 0435—1223. This single 
epoch observation, with microlensing AM^^ = —0.19 ± 0.04 and chromatic microlensing 
AMi,_R = —0.13 ± 0.08, allows us to estimate the disk size and to constrain the power-law 
index in the size- wavelength scaling of the accretion disk. We found an R-band disk size of 
= 13 ± 5 light days for a linear prior on r^, and of rf^ = 7 ± 6 light days for a logarithmic 
prior, and a value of p = 1.3 ± 0.3 consistent with the Shakura & Sunyaev (1973) thin disk 
model. Our result is in good agreement with those of other authors (Pooley et al. 2007, 
2009; Morgan et al. 2010; Blackburne et al. 2011). 

In the case of FBQ 0951+2635 and Q 0142—100, unfortunately, further observations 
are needed to completely understand the chromatic behavior of the systems, although our 
observations shed some light on their chromatic variations. The lack of chromaticity observed 
in the FBQ 0951-1-2635 NOT data is compatible with the absence of extinction, and the 
different chromatic behaviors observed at different epochs suggest that there is significant 
microlensing in the system. In the case of Q 0142—100, microlensing scenarios are also 
favored to explain the observed chromaticity, although wavelength-dependent contamination 
by the lens galaxy in the flux of image B, due to its proximity, is not ruled out. 
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SDSS J1650+4251 



Modelled image 



Residuals 



Fig. 1. — Str-y image (left) of SDSS J1650+4251, our photometric model consisting of two point 
sources (middle), and the residuals (right) after subtracting the model from the data. The peak 
residuals (white) are roughly 2% of the quasar peak. 
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Fig. 2. — Difference in magnitude vs. the inverse of the observed wavelength for the quasar images 
of the lens system SDSS J1650+4251. The filled squares correspond to our NOT data and the open 
triangles to those obtained by Morgan et al. (2003). The variation in the magnitude difference 
with wavelength is probably produced by differential extinction in the lens galaxy. An extinction 
law similar to the SMC fits the observations well, while a Milky Way extinction law (CCM) is ruled 
out. 
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Fig. 3. — Difference in magnitude vs. tlie inverse of tlie observed wavelengtli for tlie quasar 
images of the lens system SDSS 1650+4251. The filled squares correspond to the data taken at 
the NOT (2003) and the open squares correspond to data taken at the LT (2006). The trend with 
wavelength is the same in all the different epochs, although a magnitude shift due to microlensing 
may also be present. The wavelength regions corresponding to the most prominent quasar emission 
lines are indicated (from left to right: Mgll, Clll, and CIV ). The horizontal error bars correspond 
to the filter widths. 



- 22 - 



HE 0435-1223 
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Fig. 4. — Magnitude differences as a function of the inverse of the observed wavelength for 
HE 0435-1223 (NOT data on HJD 2454404). Images A and D show a clear signature of microlens- 
ing in the continuum but not in the emission lines. Moreover, the chromatic variation observed 
in image A could only be explained as a consequence of chromatic microlensing. The wavelength 
regions corresponding to the most prominent quasar emission lines are indicated (from left to right: 
Mgll, cm, CIV and Lya). The horizontal error bars correspond to the filter widths. 
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Fig. 5.— HE 0435-1223 SMARTS R-band light curve (Blackburne & Koclianel^ 2010) corrected 
for time delays. The filled squares correspond to our observation in the redshift zero Ha filter, 
which is a good match to the R-band wavelength. 
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Fig. 6. — Probability of reproducing the observed microlensing magnification and chromaticity 
of HE 0435—1223 as a function of the R-band disk size, r^, and the power-law index, p, where 
oc A^. The panels show the results for different stellar mass fractions a. The contours correspond 
to 15%, 47%, 68%, and 90% confidence intervals. Uniform (solid line) and logarithmic (dashed line) 
priors on were assumed. 
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Fig. 7. — Probability of reproducing the observed mierolensing magnification and chromaticity 
of HE 0435—1223 as a function of tlic R-band disk size, for different stellar mass fractions a, 
and normalized by the total probability in the (r^, p) grid. Uniform (solid line) and logarithmic 
(dashed line) priors on rf- were assumed. 
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Fig. 8. — Probability of reproducing the observed microlcnsing magnification and chromaticity 
of HE 0435—1223 as a function of the power-law index, p for different stellar mass fractions a, 
and normalized by the total probability in the (r^, p) grid. Uniform (solid line) and logarithmic 
(dashed line) priors on rf- were assumed. 
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Fig. 9. — Probability of reproducing the observed microlensing magnification and chromaticity 
in HE 0435—1223 as a function of the R-band disk size, r^, after marginalizing over a. The solid 

line was obtained using the prior on a from Mediavilla et al. (2009). The dashed line corresponds 
to the probability distribution without this prior. These results are for the uniform prior on r^. 
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Fig. 10. — Probability of reproducing the observed microlensing magnification and chromaticity 
in HE 0435—1223 as a function of the power-law index, p, after marginalizing over a. The solid 
line was obtained using the prior on a from Mediavilla et al. (2009). The dashed line corresponds 
to the probability distribution without this prior. These results are for the uniform prior on rf. 
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Fig. 11. — Probability of reproducing the observed microlensing magnification and chromaticity 
HE 0435—1223 as a function of tlie stellar fraction, a., marginalized over the R-band disk size 
and the power-law index p. The solid line was obtained using the prior on a from Mediavilla et 
al. (2009). The dashed Une corresponds to the probability distribution without this prior. These 
results are for the uniform prior on rf'. 
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FBQ0951+2635 




Fig. 12. — Magnitude differences as a function of the inverse of the observed wavelength for 
FBQ 0951+2635. The NOT observations (open squares) arc compared witli the CASTLES (open 
triangles) and J. A. Muhoz et al. (2011, in preparation, filled squares) HST observations. The 
wavelengths corresponding to the most prominent quasar emission lines are indicated (from left to 
right: Mgll, CIII, CIV and Lya). The horizontal error bars correspond to the filter widths. The 
solid vertical line in the Lya band corresponds to the expected position of the 2175 A feature. 
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Fig. 13. — Magnitude differences as a function of the inverse of the observed wavelength for 
Q 0142—100 at different epochs. The open squares correspond to the NOT data. The filled 
squares correspond to the data from Ehasdottir et al. (2006), and the filled triangles are the data 
from CASTLES. The wavelengths corresponding to the most prominent quasar emission lines are 
indicated (from left to right: Mgll, CIII], CIV and Lya). The horizontal error bars correspond to 
the filter widths. The solid vertical line corresponds to the expected position of the 2175 A feature. 
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Str is short for Stromgren. 

Iac#28 and Iac#29 are non-standard filters shared with the IAC-80 telescope 
(Observatorio del Teide, Tenerife, Spain). These two filters have an FWHM of 20 
nm and 9.5 nm, respectively. 
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Table 2. log of RATCam Observations 



Target Observation Date Filter Exposure (s) 



SDSS 


J1650+4251 


2006 


Jun 14 


u' 


(A = 


3500 


A) 


7 


X 


120 


SDSS 


J1650+4251 


2006 


Jun 14 


g' 


(A = 


4800 


A) 


4 


X 


100 


SDSS 


J1650+4251 


2006 


Jun 14 


r' 


(A = 


6250 


A) 


2 


X 


100 


SDSS 


J1650+4251 


2006 


Jun 14 


i' 


(A = 


7700 


A) 


2 


X 


100 


SDSS 


J1650+4251 


2006 


Jun 14 


z' 


(A = 


9100 


A) 


3 


X 


120 


SDSS 


J1650+4251 


2006 


Jun 24 


u' 


(A = 


3500 


A) 


3 


X 


120 


SDSS 


J1650+4251 


2006 


Jun 24 


g' 


(A = 


4800 


A) 


3 


X 


100 


SDSS 


J1650+4251 


2006 


Jun 24 


r' 


(A = 


6250 


A) 


2 


X 


100 


SDSS 


J1650+4251 


2006 


Jun 24 


i' 


(A = 


7700 


A) 


1 


X 


100 


SDSS 


J1650+4251 


2006 


Jun 24 


z' 


(A = 


9100 


A) 


3 


X 


120 


SDSS 


J1650+4251 


2006 


Jun 30 


u' 


(A = 


3500 


A) 


1 


X 


120 


SDSS 


J1650+4251 


2006 


Jun 30 


g' 


(A = 


4800 


A) 


2 


X 


100 


SDSS 


J1650+4251 


2006 


Jun 30 


r' 


(A = 


6250 


A) 


1 


X 


100 


SDSS 


J1650+4251 


2006 


Jun 30 


i' 


(A = 


7700 


A) 


5 


X 


100 


SDSS 


J1650+4251 


2006 Jul 7 


u' 


(A = 


3500 


A) 


3 


X 


120 


SDSS 


J1650+4251 


2006 Jul 7 


g' 


(A = 


4800 


A) 


5 


X 


100 


SDSS 


J1650+4251 


2006 Jul 7 


r' 


(A = 


6250 


A) 


4 


X 


100 


SDSS 


J1650+4251 


2006 Jul 7 


i' 


(A = 


7700 


A) 


4 


X 


100 


SDSS 


J1650+4251 


2006 Jul 7 


z' 


(A = 


9100 


A) 


4 


X 


120 



Table 3. SDSS 1650+4251 LT PHOTOMETRY 



Filter tub — ruA 



u'(A = 


3500 A) 


2.12±0.05 


g'(A = 


4800 A) 


1.90±0.05 


r'(A = 


6250 A) 


1.76±0.05 


i'(A = 


7700 A) 


1.70±0.05 


z'(A = 


9100 A) 


1.59±0.06 



Table 4. not-photometry 



Target 


Am 


Str-u 


Str-v 


Str-b 




Str-y 


Iac#28 


Ha 


Iac#29 




I-band 


SDSS 1650+4251 


iriB 


- iriA 


2.37 ± 0.09 


2.25 + 0.09 


2.12 + 0. 


,09 


2.00 + 0.09 


1.95 + 0.09 


1.90 + 0.09 


1.85 + 0. 


,09 


1.93 + 0.09 


HE 0435-1223* 


rriB 


- rriA 




0.65 + 0.02 






0.71 + 0.01 




0.66 + 0.02 






0.61+0.05 




ms 


- mc 




-0.01 + 0.03 






0.01 + 0.01 




0.04 + 0.03 






0.07 + 0.05 




TUB 


- mo 




0.15 + 0.03 






0.16 + 0.01 




0.22 + 0.02 






0.22 + 0.04 


HE 0435-1223** 


rriB 


- ruA 


0.87 + 0.14 


0.69 + 0.02 


0.81 ± 0, 


,08 


0.72 + 0.09 


0.72 + 0.10 


0.68 + 0.02 


0.67 + 0, 


,07 


0.61 + 0.04 




TUB 


- mc 


0.06 + 0.14 


-0.03 + 0.03 


0.03 + 0. 


,09 


0.01 + 0.11 


0.05 + 0.09 


0.06 + 0.02 


0.04 + 0, 


,07 


0.07 + 0.03 




TUB 


- mo 


-0.06 + 0.13 


-0.24 + 0.02 


-0.13 + 0. 


,06 


-0.21 + 0.07 


-0.20 + 0.07 


-0.23 + 0.02 


-0.24 + 0, 


,06 


-0.19 + 0.05 


FBQ 0951+2635 


TUB 


- ™-A 


1.30 + 0.05 


1.42 + 0.05 


1.37 + 0. 


,05 


1.44 + 0.05 


1.33 + 0.05 




1.36 + 0, 


,05 


1.43 + 0.05 


Q 0142-lOOt 


TUB 


- mA 


2.23 + 0.05 


2.31 + 0.05 


2.37 + 0. 


,05 


2.29 + 0.05 




2.21 + 0.05 






2.39 + 0.05 


Q 0142-100 tt 


TTlB 


- ruA 


2.24 + 0.05 




2.36 + 0. 


,05 




2.18 + 0.05 




2.37 + 0, 


,05 


2.36 + 0.05 



* Observation date: 2007 October 15. 
**Observation date: 2007 October 31. 
t Observation date: 2006 September 23. 
tt Observation date: 2006 November 17. 
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Table 5. SDSS J1650+4251 Component Positions. 



COMPONENT 


WIYN* 


NOT 






AR.A. 


ADecl. 


AR.A. 


ADecl. 


Image A 


EE 


= 


= 


= 


Image B 


0.223 ±0.002 


1.163 ±0.001 


0.24 ±0.02 


1.17 ±0.02 


Lens galaxy G 


0.017 ±0.032 


-0.872 ± 0.026 


0.015 ±0.001 - 


-0.88 ±0.07 



* Relative positions to image A obtained by Morgan et al. (2003) with the WIYN 
3.5m telescope at the Kitt Peak National Observatory. 



Table 6. HE 0435-1223 R-band Disk Size and Power-law Index Estimations for a < 0.15 





Linear Prior 


Logarithmic Prior 


a 


(light days) 


V 


(light days) 


V 


0.01 


9±2 


1.3 ±0.3 


9±2 


> 2 


0.05 


13 ±4 


1.3 ±0.3 


7±6 


1.7 ±0.3 


0.1 


7±6 


1.3 ±0.3 


7±4 


1.3 ±0.3 


0.15 


11 ±6 


1.7 ±0.3 


9±4 


1.7 ±0.3 



